Energy Reception for Mobile Radio 

By E. N. GILBERT 

(Manuscript received July 2, 1965) 

Statistical properties are derived for mathematical models of the multipath 
fading encountered in mobile radio. These properties are used to compare 
some receiving systems which use several antennas to combat fading. Par- 
ticular attention is given to a system of J. R. Pierce which has electric 
and magnetic dipole antennas and computes the electromagnetic energy 
density at a point. The statistical properties considered here include energy 
density distribution functions, correlation coefficients, and the power spec- 
trum of the energy density observed at a moving point. 

I. INTRODUCTION 

A radio signal may reach a receiver via several paths because of 
reflections from nearby objects. If the receiver is a mobile radio instal- 
lation, the received field strength may fluctuate wildly because the 
reflected waves add with changing relative phases as the receiver moves. 
J. R. Pierce* has suggested a way to combat these fluctuations by 
using three antennas. 

As background for Pierce's idea consider the standing wave pattern 
produced when a plane wave is reflected at normal incidence from a 
large wall. An electric dipole antenna moving toward the wall finds 
nulls in the electric field repeated at half -wavelength intervals. However, 
these nulls occur at maxima of the magnetic field. In fact, the total 
electromagnetic energy density ^(e|£'| 2 -f-Ai|#| 2 )is constant through- 
out the pattern. 

In Pierce's scheme, the transmitter radiates a vertically polarized 
wave. The receiver carriers a vertical electric dipole antenna and also 
a pair of loop antennas with axes perpendicular to each other and to 
the dipole. These three antennas receive the three nonzero field com- 
ponents E z , H x , and H y . The three antenna signals enter separate 
square-law detectors and the three detector outputs are added to 
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obtain } (« | E. | 2 + n I H x | 2 + n \ H y | 2 ) = rp T , the total energy den- 
sity. If the signal is amplitude modulated, the receiver may compute 
\f/ T ^ in order to achieve linear detection. 

The output of this energy receiver remains constant as the receiver 
moves through the above-mentioned standing wave pattern near a 
wall. In more complicated interference patterns the total energy density 
does fluctuate, although hopefully not as much as the electric energy 
density alone. This paper examines some superpositions of vertically 
polarized plane waves in order to compare the energy receiver with a 
receiver which observes only the electric field. Two other receivers 
are examined briefly in Section V. One is a diversity receiver which has 
two or more electric dipole antennas and a switching system to select 
the antenna with the strongest signal. The second squares and adds 
the outputs of several electric dipoles. 

Most of the analysis in this paper applies slightly more generally 
to a weighted energy detector which combines the electric and magnetic 
energy densities with weight factors 2d and 26 to obtain de\ E \ + 
bn | H | 2 (this is the energy density ^ T if d = b = |). Sections III and 
IV show that certain unequal weights have some slight advantages. 

In order to imitate the haphazardness of real mobile radio inter- 
ference patterns most of the field models which follow assume waves 
with randomly chosen amplitudes, phases, and directions of propaga- 
tion. Section IV finds probability distributions for the weighted energy 
density. At wavelengths longer than about 0.2 meters, such distributions 
might be used to predict the fraction of time that the signal will fade 
beyond the range of the receiver's AVC action. At shorter wavelengths, 
a fast automobile encounters fluctuations which have appreciable 
components at audible frequencies. Then questions about spectra 
(Section VIII) and correlations (Section VII) arise. 

Fig. 1 shows the energy density as a function of position when four 
waves superimpose. The four waves had equal amplitudes but the phases 
and directions were chosen to typify some of the random models which 
follow. The propagation directions made angles of 0°, 60°, 140°, and 260° 
which were measured clockwise away from a horizontal direction. 
Table I gives the code for interpreting the printed symbols as energy 
densities in db above the mean level. The square in Fig. 1 is 3.6 wave- 
lengths on a side. Fig. 2 uses the same four waves and the code of 
Table I to depict the electric energy density alone. It is immediately 
apparent that Fig. 2 represents a more violent function than Fig. 1. 
The peaks are higher (usually above 4 db), the valleys are deeper (often 
below —13 db), and intermediate levels are relatively scarce. 
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Fig. 1 — Total energy density of four superimposed waves. 



Fig. 2 — Electric energy density of four superimposed waves. 
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Table I — Interpretation of Figs. 1 and 2 



Symbol 


Interpretation 


blank 


below -13 db 





between —13 db and —10 db 




_10 » « -7 " 


1 


_7 » " _5 " 




-5 " " -3 ' 


2 


-3 " " -1 " 




_1 « « o " 


3 


" " 1 " 


/ 


1 " " 2 " 


4 


2 " " 4 " 


+ 


above 4 db 



II. NOTATION 

Fields will be functions of Cartesian coordinates (x,y) in a horizontal 
plane. The propagation direction of a vertically polarized wave will be 
specified by a unit vector u = (u x , u y ). Let P be the radius vector to 
a point in the (x,y) plane. At P, the following are the nonzero field 
components of a wave propagating in direction u: 

= e~ h A exp { —ifiu-P] 

~*u«A exp {— ifiU'P) 



H x = ,u 



H y = — m h u x A exp {— i0u-P] 

where 27r//3 is the wavelength and A is a complex amplitude. All fields 
depend on time through a complex factor exp iwt which will be sup- 
pressed. The factors containing the dielectric constant e and permeability 
H were inserted to simplify the expressions for energy density. When 
waves from directions u,v,w, • • • are added, their amplitudes will be 
called A (u),A (v),A (w), ■ ■ ■ . Most of this paper is concerned with the 
weighted energy density 

f(P) = d | J^A (u) exp - ipu-P | 2 



+ b | y^.A (u)u u exp — i@u-P 
+ b | £).A (u)u x exp — ifiu-P 



(1) 



The coefficients d and b in (1 ) will always be nonnegative and will satisfy 
d + b — 1. Let ypg{P) and \p H (P) denote the electric and magnetic 
energy densities at P; then \p(P) = 2dip E {P) + 2b\J/ H (P). The choice 
d = 6 = 5 makes ^(jP) the total energy density yp T (P)',^(P) can also 
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become 2\j/ E (P) or 2^ W (P) if one adopts the extreme values d = 1 or 
6=1. Another form of (1) is 

+ (P) = £4(«)A*(»)(d + bu-v) expip(v - u)-P (2) 

U,( 

where * denotes complex conjugate. 

The average level i^o about which \p(P) fluctuates, may be defined as 
the limit as R — > o° of the average of ip(P) over a circle of radius R. 
In the limit, terms of (2) with v ^ m contribute zero to the average. 
Thus, if no two propagation directions are the same the average level is 

* -HUGO p. (3) 

u 

When d = 1 or b = 1, (3) shows that the average electric and magnetic 
energy densities are each 51^0 . 

According to (3), \f/ does not depend on d and b. The influence of d 
and b on some properties of the multipath fluctuations will be examined 
in subsequent sections, and (3) guarantees that the average detector 
output remains constant as d and b vary. When making such compari- 
sons it must be recognized that the random noise received by the system 
may depend on d and b. For example, if the three antennas receive un- 
correlated noises of equal power the noise output of the detector is 
proportional to d + 26. 

III. NULLS 

Since three complex equations E M = 0, H x = 0, H y = must hold 
simultaneously at a point of zero energy density, it is not obvious when 
such a zero is possible. This section gives some examples of zeros. 

When fewer than four waves add, no two having the same direction 
of propagation, no point can be a point of zero energy density. A proof 
of this fact is given in Appendix A. However, it is not necessarily de- 
sirable to have only a small number of waves. For example, consider 
two waves propagating in directions which differ by an angle t?, say u = 
(cos \d, sin %&) and v = (cos %d, —sin ^1?). Suppose the amplitudes are 
equal in magnitude but differ in phase by 8, say 

A(u) = exp {i{ip + 5)), A(v) = exp (*p). 

At a point P = (.r, //), (1 ) and (3) show that 

*(P) = <M1 + (d + b cos &) cos (5 - 2jfy sin it?)}. (4) 

Note that ^(P) attains its minimum value ^ |1 — I d 4- b cos t? | } 



1784 THE BELL SYSTEM TECHNICAL JOURNAL, OCTOBER 1965 

along the family of lines 

(2/3 sin \&)y = 5 + multiple of *■. (5) 

The "multiple" in (5) must be odd if d + b cos # > and even if d + 
b cos t? < 0. The (positive) minimum value can be arbitrarily small if 
# is small enough. 

It is interesting to compare the weighted energy density (4) with the 
electric energy density ^g(P). When d = 1, (4) becomes 

iM-P) = !*o{l + cos (8 - 20y sin &)} (6) 

a function which vanishes along the lines (5). 

Equation (4) may be used to help decide a good choice of the coeffi- 
cients d,b for a detector. Imagine the two waves produced at random in 
such a way that the angle t? and relative phase 8 are independent random 
variables, both having probability density (2ir)~ l in the interval (0, 2t). 
One wants \f/(P) to fluctuate as little as possible. The variance of the 
random variable \f/(P) is one measure of fluctuation. From (4) one ob- 
tains a variance 

E{\f(P) -h\ 2 } -**"{* + W 

which has its minimum when d = -|, b = § . Alternatively, one might 
prefer to pick d and b so that the expectation of the minimum value 
1 — | d + b cos t? | of ^(P) is as large as possible. This condition re- 
quires d and b to minimize E{ \ d + b cos # | ) . The calculation given in 
Appendix B shows that the minimizing d and b are d = 0.40, b = 0.60. 
Although both criteria suggest that the magnetic energy density be 
weighted more than the electric energy density, both minima are so 
broad that an energy detector with d = b = % does almost as well. 

It is possible to have \f/ (0 ) = when waves from four different direc- 
tions superimpose. For example, take the propagation directions along 
the ±z, zky coordinate axes 

u = (1,0), v = (0,1), -u = (-1,0), -v = (0, -1) 

and let the amplitudes be 

A(u) = A(-u) = 1, A(v) = A(-v) = -1. 

At the point P = (x, y), 

4>(P) = \fi {d(cosPx - cos/3*/) 2 + 6 (sin 2 #r + sin 2 @y)\. 

Zeros are spaced \/y/2 apart in this pattern if b ^ 0. However the 
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electric density is much worse, 

MP) = W«w/te - cos /ft/} 2 . 

The zeros of \f/n(P) are not isolated, they occupy two orthogonal families 
of parallel lines. 

IV. ENERGY DISTRIBUTION FUNCTIONS 

Three real parameters specify a wave, say the angle t? between u and 
the x axis, and the modulus | A (m) | and phase of the complex amplitude 
A(u). This section discusses some models which pick at random the 
3A7 parameters of N interfering waves. In every case the TV phases are 
chosen independently with constant probability density (27r) _1 in the 
range (0,2t). As a result, the models are stationary with respect to 
translations of the (x,y) coordinate system. In particular the probability 
distribution function of \p(P) is the same for all points P; to simplify 
the analysis take P = 0, the origin. The distribution function 

F(*) = Prob{*(0) ^ M 

is the probability that the detector of a mobile radio station produces 
an output less than \p. In this section, the N waves have roughly the 
same statistical properties. By contrast Section VI considers a model in 
which one of the waves represents a strong wave direct from the trans- 
mitter. 

In the first model the number of waves is N ^ 3. The N propagation 
vectors U\ , • • • , u N are not random. They are equally spaced around the 
unit circle; u k makes angle 2-irk/N with the positive .r-axis. Each complex 
amplitude A (u) will have the form -4 (u) = R(u) + il (u) where the 
2N real numbers R(it\), • ■ ■ > R(un), I(Mi)t ' " > I(i*n) are supposed 
independent Gaussian random numbers with mean and variance 1. 
Another way to obtain the same random process is to pick moduli 
| ^4 (wi) | , • • • , | A (u N ) | independently from a Rayleigh distribution 
and the N phases independently with constant density (2ir) _1 in the 
range to 2ir. 

According to (3), the average of \p(P) over the plane is 

*o = XX» + £/ 2 (t0 

U II 

so that the expected average is 

= Bfa) = 2N. (7) 
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Appendix C derives the distribution function 

F<&) - 1 -c 2 exp (-*'/d) 

+ (c - 1) (c + 1 + 2^/d} exp (-2*'/6) 

where c = 2d/ (2d — b) and ^ = ^/$. Note that the number of waves 
N enters (8 ) only through the normalizing factor $ = 2N. The distribu- 
tion function for the total energy is (8) with c = 2. From the limiting 
cases d = 1 and d = of (8 ) one obtains distribution functions for the 
electric and magnetic energy densities. 

Prob {^(0) £ {WW} = 1 - exp (-*') 

Prob {^(0) ^ (§#)*'} = 1 - (1 + 2*') exp (-2^). 

Curves ^JV and 7W in Fig. 3 show the distributions of electric and total 
energy densities plotted in db above their respective mean average 
levels. The electric energy density has a much higher probability of 
being small than the total energy density. A curve for the magnetic 
energy density will be given in Section V. 

The distribution obtained for \f/ E (0) implies that the electric field 
strength | E z | has a Rayleigh distribution. In this respect, the model 
agrees with some experimental data of W. R. Young (see in particular 
his Fig. 5). 

For small values of \f/, (8 ) becomes 

f&) = dim'H ••• 

with missing terms • • • of order 0(^' 4 ) (see (24)). To make small values 
of ^ as unlikely as possible, one may minimize d~ l b~ 2 by picking d = \, 
b = f . Recall that these values had another minimizing property in 
Section III. Again the advantage over using d = 6 = \ is slight. If the 
curve for d = $, b = f were plotted in Fig. 3 it would lie about J db to 
the right of the total energy density distribution curve. Equation (8) 
becomes indeterminate when d = f, b = § . However, in this special case, 
i^(0) has a chi-squared distribution of six degrees of freedom 

Fty) = 1 - (1 + 3*' + |(3^') 2 ) exp - 3*'. 

Part of the variability of ^(0) comes from the randomness of the 
average value \f/ . For any particular choice of the wave amplitudes, 
^o will not be exactly #; the distribution of ^(0)/^o might have been 
more relevant. However, \f/ has a chi-squared distribution with 2N de- 
grees of freedom and so has high probability of being close to $, say 
within 0.5 db, especially when N is large. 
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ENERGY DENSITY (IN DECIBELS ABOVE MEAN VALUE) 

Fig. 3 — Probability distribution functions for energy densities. 

The simple form of the distribution (8) results from the special ran- 
dom process which picks the amplitudes and directions. One might 
prefer to choose directions independently at random with probability 
d/2ir of making an angle less than 1? with the z-axis. Other amplitude 
distributions also suggest themselves. It seems reasonable to continue to 
insist on independent amplitudes with random phases but one might 
use equal moduli | .4 (w) | or another modulus distribution instead of the 
Rayleigh distribution. Undoubtedly Fty) will be a more complicated 
function of N in these cases. However, (8) must still apply in the limit 
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of large N. For example, let M be an integer approximately equal to 
N . Let Vi , • ■ • , v M be M unit vectors equally spaced around the unit 
circle. For each wave direction u find the vector u in the list Vi , V2 , • • • , 
v M which approximates u as closely as possible (and so to within angle 
ir/M). If N is large, one makes only a small error in i^(0) by replacing 
each true direction u in (1) by its approximation u . This approxima- 
tion replaces the waves with random directions by waves with M equally 
spaced directions. For i = 1, • ■ • , M the approximating waves with 
direction y, add up to a single wave, say with amplitude A (»<). The 
central limit theorem shows that the real and imaginary parts of A (»<) 
have approximately Gaussian distributions when N is large. Then the 
assumptions leading to (8) hold again. 

Limiting results may be misleading if applied when the number of 
waves is small, as may be typical in mobile radio. Curves T2 and E2 of 
Fig. 3 show distributions for the total energy density and electric energy 
density (plotted in db above their mean values) for a superposition of 
two waves with equal moduli, random phases, and random directions. 
The distribution of total energy density yJ/ T (0) was obtained numerically 
using (4) with d = b = £; ^ r (0) was evaluated for 200 equally spaced 
values of t? and 200 equally spaced values of 8. A histogram of the 40,000 
numbers was compiled to get the distribution. The electric energy dis- 
tribution is easily derived from (6): 

Prob [$b(0) ^ (hfo)t'\ = 7T _1 arc cos (1 - +') 

(see Margaret Slack 4 for the electric energy distribution when other 
numbers of random waves of equal moduli combine). The curves show 
that the case of two waves of equal moduli is much worse for mobile 
radio than the case of waves with Rayleigh distributed moduli. Never- 
theless total energy detection is again much better than electric energy 
detection. 

When more than two waves of equal moduli, random phases, and 
random directions combine, the total energy density \f/ T (0) depends on 
many random parameters. To find its distribution function by a nu- 
merical integration of the kind used for two waves would be much too 
costly. A computer experiment was used instead. Using pseudo-random 
numbers to pick phases and directions, the computer generated a se- 
quence of field components for independent plane waves. After com- 
puting each new wave the computer found for N = 2, • • • , 10, the 
energy density in the sum of the N most recently computed waves. 
After 10,009 waves, the computer had compiled histograms of the energy 
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densities in sums of 2, 3, • • • , 10 waves, each based on 10,000 random 
samples. The same wave appeared in N consecutive sums of N waves; 
then samples closer together than N were not independent. However, 
the estimate of 7' 1 (if/) is at least as good as if the experiment had 10,000/iV 
independent samples. 

Table II summarizes this experiment and compares the observations 
with theoretical predictions based on the curves in Fig. 3. The numbers 
observed agree surprisingly well with (8) even when N = 3. A compari- 
son of the theoretical and observed numbers for N = 2 gives an idea of 
the accuracy of the experiment. 



V. RECEPTION USING if ELECTRIC DIPOLES 

Let m vertical dipole antennas be placed at points Pi , • • • , P m . 
Using switched diversity reception the received signal energy density is 

yp D = max (fr(Pi), ■■• i **(Pm)}> 

Another possibility is to square the antenna signals and add them 
(additive diversity) ; then the detector output is 

ia = MPi) + MP*) + • • • + MPm). 

Pi , • • • , P m will be assumed spaced so far apart that the fields at 
these points may be considered independent random variables. If there 
are N waves generated at random by the first model of Section III, then 
each term ^«(P,-) is a random variable with the chi-squared distribution 
of two degrees of freedom and mean N. 



Table II — Fraction of Sums of N Waves of Energy Density 
^ f. Sample Size is 10,000 





N = 2 


N = 2 


N = i 


N = 5 


N - 10 


large W 


in db 


theor 


obs 


obs 


obs 


obs 


Eq. (8) 


-16 


0.0057 


0.0066 


0.0002 








0.00008 


-14 


0.0091 


0.0084 


0.0002 


0.0003 





0.0003 


-12 


0.0144 


0.014 


0.0004 


0.0003 


0.001 


0.0011 


-10 


0.0232 


0.025 


0.0004 


0.002 


0.003 


0.0042 


-8 


0.0376 


0.038 


0.006 


0.009 


0.010 


0.0144 


-6 


0.0614 


0.063 


0.030 


0.033 


0.047 


0.0459 


-4 


0.1037 


0.104 


0.105 


0.120 


0.117 


0.1301 


-2 


0.1851 


0.185 


0.227 


0.270 


0.296 


0.3103 





0.5000 


0.500 


0.530 


0.574 


0.584 


0.5869 


2 


0.8908 


0.891 


0.882 


0.853 


0.849 


0.8484 


4 


1 


1 


0.995 


0.989 


0.979 


0.9742 


6 


1 


1 


1 


1 


1 


0.9986 
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The distribution function for \f/ D is 

F D (f) = Prob (& D ^ *) = II Prob {*,(P<) ^ *} 
= {1 - exp (-HN)r- 
The mean of ^ D is 

f - jT {1 - 7,(*) }<** = g (j) (-l) w /fe 

For m = 1, 2, • ■ • , 5, $ D is N, 3AT/2, lliV/6, 25N/12, 137W/60. Then 

F D (f) = {1 - exp (-3^/2i? D )j 2 whenm = 2 

FdM = jl - exp (-ll^/6i? D )) 3 whenm = 3, etc. 

\f/ a has the chi-squared distribution with 2m degrees of freedom and 
mean Nm. Fig. 4 shows the distribution function of ^ s with m = 2, 3, 
5, and 8. The curves for the distributions of ty D and \l/ s lie very close and 
so the curves for \p D were not added to Fig. 4. 

When m = 3 the distribution function of \{/s is exactly the same as 
the one for the weighted energy density ^(0) with d = \ and b = f. 
As noted in Section IV, this distribution function is slightly better than 
the one for the energy density (d = b = %) at small values of ^. 

In Section IV, the two magnetic field components H x (0) and H v were 
found to be uncorrelated. Then the distribution function of the magnetic 
energy density at zero follows the curve labeled m = 2 in Fig. 4. 

VI. STRONG DIRECT WAVE 

Section IV presented extreme cases in which the waves are all roughly 
of comparable strength. In this section another wave, stronger than the 
others, will be added to represent a "direct" wave of amplitude R. It is 
no longer easy to derive P (^ ) exactly. However, the asymptotic form of 
F (^ ) for small values of \f/ is derived in Appendix C using the first model 
of Section IV. This result (25) assumes a convenient form in terms of 
the quantity 

a- = 2N/$ = 1 - R 2 /$ 

which represents the fraction of the expected weighted energy density $ 
contributed by the N scattered waves. When N and R are expressed in 
terms of $ and a the final result is 

Fty) = (D<TVV 8 (*/*) 8 exp - 3(<r' - 1) (9) 

approximately for small \f/. 



ENERGY RECEPTION FOR MOBILE RADIO 



1791 



10-2 



10-3 



10"" 











/£?^~~ 
















































































/ 












/I 












I I 












1 1 












III 










j J 






m =2/ 


3/ 


,/ /. ELECTRIC DIPOLE 
3 / 1° ANTENNAS 








/ / 












/ / 












/ / 












/ / 








/ 


1 1 












1 










^ 










1 


1 










J 


J 










/ 


1 











10-5 

-20 -15 -10 -5 5 10 

Fig. 4 — Probability distribution function for the sum yp B of the electric energy 
densities at m points. 

The case R = 0, or <r = 1, was discussed in Section IV. When < 
a < 1, (9) contains the extra factor 

a~ exp — 3(<r -1 — 1) 

which is less than 1 and approaches with decreasing a. One concludes 
that F (^) is then smaller than the value given by (8) ; i.e., deep minima 
tend to be less frequent when a direct wave is present. For example, if 
scattered waves account for only half the received weighted energy 
density (a = §), the probability of receiving less weighted energy 
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density than \p is only 8 exp — 3 = 0.40 times the probability (8) for 
a = 1. 

VII. CORRELATION COEFFICIENTS 

This section finds correlation coefficients between various pairs of 
energy densities. The correlation coefficient between ^(0) and yp(P) 
indicates whether a receiver traveling from to P will find very different 
weighted energy densities at the two points. The correlation coefficient 
between ^b(O) and rf/siP) might be used to decide whether and P are 
good locations for two electric dipole antennas in a diversity system; 
one would want low or negative correlation. For similar reasons, corre- 
lation coefficients involving the magnetic energy density ^«(P) = 
\I/(P) — ^b(P) are interesting. 

The waves in this section are produced by a random process slightly 
more general than the one used to get Table I. The N propagation 
directions and N phases are chosen at random and independently as in 
Section IV. The moduli | A (wi( | , ■ ■ • ,\ A (u N ) | are now chosen in- 
dependently from a common probability distribution. It is not necessary 
to know the distribution in detail. Only the expected values of the second 
and fourth powers E(\ A \ 2 ), E(\ A \ ) enter into the correlation co- 
efficients. To facilitate comparisons with Section IV, take E( \ A \ 2 ) = 2. 
Then the expected average weighted energy density is again 

$ = 2N. 

R. H. Clarke has also used this random process in an unpublished study 
of some different correlations. 

It will be convenient to express the fourth power moment as 

E{\A\*) = 4+2 2 , 

so that 2 2 is the common variance of the squared moduli. When the 
moduli are all the same, | A (w,) \ = 2\i = 1, • ■ • , N, and the variance 
S 2 is zero. When moduli have the Rayleigh distribution, S 2 = 4. 

All the correlation coefficients of interest will be obtained as special 
cases of a single result. Consider two weighted energy densities. 

in = 2# fi (0) + 2ty ff (0) 

and 

f 2 = 2D+ E (P) + 2B^ H (P). 
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Appendix D proves 

Efah) - E(\l>i)Ety 2 ) = N2 2 + 4N(N - 1) {dDJ o 2 (0r) 

+ (dB + bD)J 1 2 (fir) (10) 

+ lbB(jf(pr) + J?(fir))\ 

where r = | P | . When D = d, B = b, and P = 0, then ^ 2 = ^i and 
(10) becomes the variance of i^i 



Var ft = NZ 2 + (4d 2 + 2b 2 )N(N - 1). (11) 

Likewise, 

Var^ 2 = N2 2 + (4Z> 2 + 2B 2 )N(N - 1). (12) 

The coefficient of correlation between \f/\ and ^ 2 is 

P = [EbHd - Edh)E(rh)}/ {Var fc Var fc} *, (13) 

which may be evaluated using (10), (11), and (12). The case of equal 
moduli (S 2 = 0) is especially simple because then (assuming N > 1) 
the factors N(N — 1 ) cancel out and p does not depend on N. This is 
the only case in which p — > as r — > <x> ; when S 2 > the average 
weighted energy density \f/ is uncertain and so the energy densities re- 
main slightly correlated even at points far apart. When N is large this 
residual correlation is small and p approaches its value for equal moduli. 
By choosing special values for d, b, D, and B, one can obtain correla- 
tion coefficients of special interest: 

Ptt = !3J« 2 (/3r) 4- 4 I / 1 2 (0r) 4- J 2 2 (fir)\/S 

Pbb = Jo 2 (fir) 

Pbb = Jo 2 (fir) + Jiifir) 

Pkt= (DMtffor) + Jx(fir)} 

PEH = 2' J {(fir) 

Pht = 3~*Uo 2 (/3r) + 2J{(fir) + J?(fir)\. 

Here the subscripts E,H,T indicate the kind of energy, electric, mag- 
netic, or total, at the two points. For the sake of simplicity the coeffi- 
cients have been given only in the special case of equal moduli (S 2 = 0) 
or in the limit of large N. Some of the more interesting coefficients are 
plotted in Fig. 5. 
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1.0 




\/2 X 3/2 X 

DISTANCE, T, BETWEEN POINTS 

Fig. 5 — Coefficient of correlation between the energy densities at two points 
separated by distance r. 



VIII. SPECTRA 

When a receiver moves with constant velocity vector V the received 
energy density \f/(V t) is a random function of time. If one assumes the 
model of Section VII, the autocorrelation function of \p(Vot) is known 
and hence its power spectrum may be found. The power spectrum of 
^(FoO gives some idea of the frequencies at which the fluctuation noise 
is likely to be strong. In particular cases, depending on the way that 
the modulating signal is to be extracted from ^(P), the spectra of other 
functions may be more important. For example, if an AM system is 
used, one might prefer to know the power spectrum of \f/ (Vot). For 
another kind of fading spectrum see J. F. Ossanna . Ossanna combines 
two random waves and derives the spectrum obtained at the output of 
an envelope detector receiving the electric field. 

For purposes of comparing power spectra it is convenient to nor- 
malize them to make the total power unity. Appendix E takes the 
Fourier transform of E{\l/(0)\l/(V t)}/E {^(O)} to obtain a normalized 
spectrum. In order to keep formulas simple, Appendix E and this sec- 
tion consider only the case of large N. 

The normalized power spectrum of \f/(V t) contains a spectral line at 
zero frequency which represents the carrier or desired signal. The power 
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179, f 



in this line is 1/(1 + r/ 2 + %b 2 ). Again the choice d = %, b = f maxi- 
mizes this power. The rest of the spectrum is fluctuation noise distributed 
with a spectral density function s(f). Fig. 6 shows this spectrum for the 
total energy density and for the electric energy density. Both spectra 
vanish when / is larger than a cutoff frequency 

/o = 2 | F |/X. 



(14) 

Note that / is the frequency of the fluctuations in electric energy den- 
sity observed by a vehicle moving toward the wall in the interference 
pattern described in Section I. When < / ^ /o , the spectral density 
has an analytic expression (32 ) in terms of the complete elliptic integrals 
K(x) and E(x). Let v = ///„ . Then 



s(f) = 



16 



(3 - v 2 )K((l - v 2 ) 1 ) - 2(2 - V 2 )E((\ - v 2 ) h )} (15) 



33tt 2 /o 
for the total energy density and 

a(f) = K{(1 - v 2 )*\/(*%) 



(16) 



for the electric energy density. 

The fluctuation noise \J/ T (Vot) appears to be less troublesome than 
\f/x(V( ,t) both because it contains less total power away from the carrier 
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Fig. 6 — Power spectra of the energy densities i/<(VoO and if>t;(V t) observed 
by a vehicle moving with constant velocity Vo . 
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and also because its spectrum is concentrated more toward low fre- 
quencies. As Fig. 6 shows, the spectrum of \l/ T (Vut) goes to smoothly 
at/ = /o while the spectrum of \f/ B (Vot) remains at a high level until it 
drops to discontinuously at / = /o . 
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APPENDIX A 

The impossibility of producing a zero by adding fewer than four waves 

In what follows, waves must have nonzero amplitude and no two 
waves may have the same propagation direction. It is clearly possible 
for two or three waves with same direction to cancel if their amplitudes 
add up to zero. 

The condition for a zero at the origin (P = 0) is 

o = *(o> = 1 5> («) r + 1 5> ( w K r + 1 1> (^k r 

u u u 

or simply 

= 2>(w) = J^A(u)u. (17) 

Consider first the case of two waves in different directions u and v. Then, 
(17) becomes 

A(u) +A(v) = 0. 

A (u)u -\- A (v)v = 0. 

The second (vector) equation requires that the unit vectors u,v be 
colinear. Since v cannot equal u, v = — u. Then 

A (u) + A (v) = 

A{u) - A (v) = 0, 

a system with no solution except the trivial one A (u) = A(v) = 0. 

When there are three waves with directions u,v,w, one may eliminate 
A (w) from the system (17) to get 

A(u)(u — w) + A(v)(v — w) = 0. 
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A (u) cannot be zero. Then 

u = w + a(v — w) 
where a = —A (v)/A (u). Since | u | z = | w \ = 1, one finds 

= 2aw (v - w) + a 2 \ v - w | 2 . (18) 

Also, 

= 2w (v - w) + | v - w | 2 (19) 

follows similarly from v = w + (v — w) and | v \ 2 = | w | 2 = 1. Use (19) 
eliminate 2w- (v — w) from (18) and get 

| v — w \ 2 a(a — 1) = 0. 

Since v j± w and since fl^O (otherwise A(v) = 0), a = 1. However, 
ifa= 1, u = w -\- l(v — w) = v, a contradiction. 

APPENDIX B 

Weights which maximize the expected minimum value of\f/(P) 

In the interference pattern (4) for two random waves, rf/(P) attains a 
minimum value 

^min = '/'ofl - I d + b cos i? | }. 



The expected value of ^ min is 

E(+ min ) - to [l - fa)' 1 f I d + b cos & I dtfj . (20) 

The evaluation of the integral in (20) requires two cases. First, if b tz 
| g: d, \ d -{- b cos tf | = d + 6 cos i? and 

E(+ min ) = U, (b ^ d). (21) 

Second, if d g \ ^ 6, let tf = cos -1 (fZ/6). Then 

| d + ?> cos i? | = f d + 6 cos # when | t? ^ tt - d 

[ — d — b cos t> otherwise 
and 



E(+, ) = Mb + (2/t)(cW -V6 2 - d 2 )!, (<* ^ &)• (22) 

Fig. 7 shows E(rl/ m - in ) plotted vs d. There is a broad maximum near 
d = 0.4, b = O.G. 
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Fig. 7 — Ety ml «) and E(^ min )/(d + 2b). 

As mentioned in Section II, the received noise power will depend on 
d and b. One might prefer to maximize the expected minimum signal-to- 
noise ratio. If the three antennas have uncorrelated noises of equal 
powers one would then maximize E(\f/ min )/(d + 2b). The dashed curve 
in Fig. 7 shows that d = 0.45, b = 0.55 for the maximizing detector. 

appendix c 



The energy distributions (8) and (9) 

It is convenient to have a special notation for real and imaginary parts 
of the field components at P = 0. S will always be a real part, J will 
always be an imaginary part. Subscripts 1, 2, or 3 on S or J denote the 
field, either £* E 2 , \x h E x , or ti h H y . For example, the imaginary part of 
M*#* is 

J i = ^2u u I(u). 

u 

Each of the six components Si , Ji , S 2 , J 2 , S s , J 3 is a linear combina- 
tion of the 2N Gaussian variables fi(%), • • • , I(u N ). Then these six 
variables have a joint Gaussian distribution, which is determined en- 
tirely by its 21 -second moments. 

All second moments of the form E(SiJj) are zero. This follows be- 
cause E(R(u)U(v)) = E(R(u))E(I(v)) = for all W 2 choices of 
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u,v. Two other typical second moments are: 

E(SM = £{-£* (*)#(»)*! 

11,1 

= -I>*B0B 2 (w)) 
- - J>. 

= - 2 cos (2*k/N), 



and 



The identity 



E(StS z ) = E{-^R(u)R(v)u u v x \ 

u,« 

" „ 1 - exp (i'M) 

>. exp (t/a) = e — t^tt- 

V 1 — exp (il) 



can be used to prove that both are zero. In the first case set t = 2r/N 
and take the real part (recall that N ^ 3 is assumed). In the second 
case take t = 4jt/N and take the imaginary part. In like manner, one 
eventually finds that the only nonzero moments are 

EOS! 2 ) = Wi 2 ) =Sl-iV 

u 

E(S?) = E(J.I) = L sin' ** = L (± - ieoa 4 -^) = itf 

E(Sa ! ) = W,') = S cos* ** = S (i + i cos **) = iJV. 

Note that these formulas hold only because N ^ 3. The cases N = 1 
and 2 are different, having E(S*) = E(J 2 2 ) = and E(S 3 2 ) = E(./ 3 2 ) - 

AT. 

The six parts of the field components are independent Gaussian 

variables with joint probability density function 

HN*Y* ex P (-{Si 8 + J * + 2 < & ' + J « f + *' + J * 2 )!/ 2A ^- < 23 > 

Now note ^(0) = dtj -\- bU where *i = S* + J? and h = <S 2 2 + J 2 2 + 
5 3 2 _j_ / 3 2 ar e two independent variables with chi-squared distributions. 
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Then the desired distribution function is 
Fty) = Prob {dh + bt 2 ^ f) 

= %N~* / feexp {-h/N) / exp {-t 1 /2N)dkdt- 1 . 

An elementary integration produces the final result (8). 

The asymptotic form of F(\f/) for small \J/ can be obtained by differen- 
tiating (8) or, more simply, by the following argument. According to 
(23 ) the joint probability density of Si , Ji , *S 2 , J% , S 3 , «/ 3 at the 
origin is \ (Afar) . The inequality ^(0) ^ ^ defines a small ellipsoid 

d(S{ + J{) + 6(& 2 + J2 + S 3 2 + Ji) ^ + 

about the origin. This ellipsoid has six-dimensional volume 

(7r76)<rVV 8 . 

The probability that (Si , Ji , • ■ • , J 3 ) lies in this ellipsoid is, apart 
from terms of higher order, 

fu,) = KtfT)- i (»7d)«rvy 

fty) = (i)<^ 1 &~ 2 (*/0)^ (24) 

In Section VI, an additional wave, stronger than the others, was 
added to represent a "direct" wave. Let the direction u of the direct 
wave be along the a>axis and let its amplitude be A (u Q ) = R, a given 
real number. With Si , S2 , • • • , J3 denned again to include the random 
fields only, the weighted energy density is 

*(0) = d { (R + Si) 2 + Ji] +&{&■+ /•' + (B + S3) 2 + ^ 3 2 } 
with mean $ = i£ 2 + 2iV. The asymptotic form of the distribution 
function for 1^(0) may be derived in the same way as (24). Now the 
six-dimensional ellipsoid ^(0) ^ ^ of volume (w 3 /Q)d~ 1 b~ 2 4 /3 is centered 
on the point (—R,0,Q,0, — R,0). Equation (23) gives the probability 
density at that point and hence the result 

F(+) = f*/[12N*db 2 exp (SR 2 /2N)} (25) 

approximately for small \J/. 

APPENDIX D 

Correlation coefficients 

To prove (10) write, 

E(ifah) = 4dDEtym(0)to(P)) + 4dBEfo(0)fB(P)) 

(26) 

+ 4bDE(f B (0)MP)) + 46JW?(^(0)^(P)). 
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The four expectations on the right are: 

E\^,(0)^ K (P)\ = N(N - 1)|1 + ./o 2 (/3r)( + \NE(\ A | 4 ) (27) 

E{+A0)+ H (P) = E\+„(0)+ B (P)\ 

(28") 
= N(N - 1)|1 + y, 2 (/3r)} +\NE(\A | 4 ) 

E{f„(o)+ H (P)) = n(n - on + w^m + \j;m\ 

+ \NE{\A 

The proofs of (27), (28), (29) are alike. Only E\t E (0)+ H (P)} will be 
derived in detail. 
Begin with 

2^(0) = | £a(«)| 2 



(29) 



Wh(P) = I L A- (tOC/ exp - iflU-P | 2 . 
Then, 

MWAP) 

= ££A(ii)A>HM'.4(t/)A*(F)exp#(F - U)-P (30) 

u.v.U.V 

with the summation variables u,v,U,V ranging over all N ways of 
picking four vectors from U\ , ■ ■ • , u N . Now take the expectation of 
both sides of (30). Most of the AT* terms in the sum have zero expecta- 
tions because E(A(u)) = E(A 2 (u)) = 0, and because A(u), A(v) 
are independent when u ^ v. Nonzero expectations can come from 
terms of three types: 

(\) N(N — 1 ) terms with v = u, V = U, and U ^ u, 
(ii) N(N - 1) terms with V = u, v = U, and U ?* u, 
(iii) N terms with u = v = U = V. 

The expectation of each term of type (0 is 

E(\A(u)\ 2 U-U\A(U){ 2 ) = lE(\A(u)\>)E(\A(U)\ 9 ) 

= |X2X2-1. 

All N(N - 1) terms contribute N(N - 1) to (28). 
The expectation of each term of type (ii) is 

\E\ | A (tt)|'| A (I7)|*i*- U exp i/3(u - U)-P\ 

= E{u-U expip(u - U)-P\ 
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Now let # and 6 be the angles which u and U make with the vector P 
so that u- U = cos (d — 0) = costf cos 8 -f sint? sin 6 and (u— U)-) 
— r cos § — r cos 9. The expectation sought is 

E {cos (i? — B) exp ipr (cos # — cos 0) | 

(I f 2 ' 
= )\ I cos tf exp (ifir cos &)d&/2w 



+ 



/•2t 

/ sin d exp (ifir cos &)d§/2ir 



The two integrals may be recognized as Fourier coefficients in the well- 
known series 

00 

exp {ifir cos &) = J a (fir) + 2 £ i n J n (fir) cos wtf. 

n = l 

Then each of N (N — 1) terms of type (ii) contributes J? {fir) to (28). 
Each term of type (iii) is \ \A(u)\*; the total contribution to (28) 
of all TV terms is \NE{\ A\*). 

APPENDIX E 

Spectra 

The normalized power spectrum of \J/{Vot) is a Fourier transform of 
E(yp(Q)yp{Vd))/E{yp 2 (Q)). The expectations are obtainable from (10) 
with D = d, B = b, P = V t and from E(f(0)) = E(t(Vot)) = 2N . 
When N is large, one seeks the transform of 

1 - dW(fir) + 2 db J?(fir) + hb\j Q \fir) + J?(fir)) ,„. 
1 + d 2 + ifc 2 K6l) 

with r = | Vo | t. 

The constant term in (31) represents the spectral line described in 
Section VIII. The remaining terms may be transformed using the equa- 
tion 



r 2 \ p 

I Jn(x) cos xy dx = \ 
10, 



P n - h {\if - D, < y < 2 



2 < y < <x> 



(Ref. 2, Erdelyi, et al, p. 46, transform 21), where P„_i(m) is the Le- 
gendre function of order n — \. 
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This result is applied in the form 

f Jn(P | Vd | ) COS 2*ft dt = (-l) n P n -i(2u 2 - l)/(7T/o) 

for ^ v ^ 1 (recall (14) and v = f/f ). One then obtains an expression 
for s(f) which involves Legendre functions of orders — §, \, and |. 
This expression was not suitable for computing because there was no 
available table of Legendre functions of fractional order. However, the 
Legendre functions of half-integer order can be expressed as complete 
elliptic integrals by the following identities: 

P_,(s) = (2/*)£{(* - hx)H 

Pi(x) = (2/ir)\2E((h - ir) 5 ) - K{{\ - \x)*)\ 

(m + l)P m+1 (x) = (2m + l)xP m (x) - mP m ^(x). 

(Ref. 1, Abramowitz and Stegun, Eqs. (8.13.1), (8.13.8), (8.13.11)). 
When the Legendre functions are replaced by elliptic integrals the Bessel 
function terms of (31 ) transform into 

2{(3-4i/ ! / 2 )K((l- , 2 )*) 

+ [(8^ - 4)6 2 - 12 db]E Ul - " 2 )*)} (32) 

s(/) - ~" 3tt 2 /o(i + ^ + m 

The results (15) and (16) are special cases of (32). 
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